Interaction Corrections to Two-Dimensional Hole Transport in Large r s Limit 



(N 
O 

o 

(N 

o 

Q 



i 

C 

o 
o 



(N 
> 

in 

O 
<N 
O 

-i— > 
c3 



i 

C 

O 
o 



X 



H. Noh, 1 M. P. Lilly, 2 D. C. Tsui, 1 J. A. Simmons, 2 E. H. Hwang, 3 S. Das Sarma, 3 L. N. Pfeiffer, 4 and K. W. West 4 

1 Department of Electrical Engineering, Princeton University, Princeton, New Jersey 08544 
2 Sandia National Laboratories, Albuquerque, New Mexico 87185 
Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park, Maryland 20742 
' Bell Labs, Lucent Technologies, Murray Hill, New Jersey 07974 
(Dated: February 6, 2008) 

The metallic conductivity of dilute two-dimensional holes in a GaAs HIGFET (Heterojunction 
Insulated-Gate Field-Effect Transistor) with extremely high mobility and large r s is found to have 
a linear dependence on temperature, consistent with the theory of interaction corrections in the 
ballistic regime. Phonon scattering contributions are negligible in the temperature range of our 
interest, allowing comparison between our measured data and theory without any phonon subtrac- 
tion. The magnitude of the Fermi liquid interaction parameter Fg determined from the experiment, 
however, decreases with increasing r s for r a > 22, a behavior unexpected from existing theoretical 
calculations valid for small r s . 

PACS numbers: 73.40.-c,71.30.+h,73.40.Kp 



In two-dimensional (2D) charge carrier systems, it is 
well known that any amount of disorder in the absence of 
interactions between the carriers will localize the carri- 
ers, leading to an insulator with zero conductivity (a) as 
the temperature (T) is decreased to zeroEl. Recent exper- 
iments on high mobility dilute 2D systems, on the other 
hand, have shown a "metallic" behavior at low T, char- 
acterized by an increasing a with decreasing T, and an 
apparent metal-insulator transition (MIT) as the carrier 
density is loweredo. There are three important energy 
scales in these systems. The first two are the Fermi en- 
ergy and the interaction energy. Their ratio, which is r s , 
is around 10 or higher for the systems where the MIT 
is observed, implying that interaction must be playing a 
role. The other energy scale is related to the disorder in 
the system given by K/t, where r is the elastic scattering 
time. It has been found from more recent experiments 
that disorder is also playing a significant role. In par- 
ticular, the critical density [n c for electrons and p c for 
holes), above which a system shows the metallic behav- 
ior, is foundr-jto decrease when disorder in the 2D system 
is decreased!! 

An important question is whether this apparent metal- 
lic state is truly a new ground state of the 2D charge 
carriers or simply a novel finite temperature behavior of 
the 2D gas, since all experiments are done at finite T. 
What is measured in such experiment is the temperature 
coefficient, da/dT. The metallic behavior evinced by the 
observation of negative da/dT at finite T does not nec- 
essarily mean, however, a true metal with nonzero con- 
ductivity at T = 0. Recently, Zala et ala calculated the 
Fermi liquid interaction corrections to the conductivity 
in the asymptotic low temperature regime {T /Tp <C 1 
where Tp is the Fermi temperature), and pointed out 
that the metallic behavior seen in the high mobility sam- 
ples could be understood by taking into account of in- 
teraction corrections in the "high temperature" ballistic 
regime (ksT 3> H/t). They found that the conductivity 



of interacting 2D carriers changes linearly with T in the 
ballistic regime, Tp 3> T ^> H/keT, with the sign as well 
as the magnitude of da/dT depending on the strength 
of the interaction, while in the low temperature diffu- 
sive regime {ksT <C h/r) the conventional logarithmi- 
cally changing conductivityO is recovered. A linear de- 
pendencft-.of a on T has also been predicted in earlier 
theoriesoQ based on temperature dependent screening, 
and this screening contribution is included in the theory 
by Zala et al.. 

Experimentally, however, it is not straightforward to 
identify the interaction corrections unequivocally in the 
ballistic regime for two main reasons. First, scattering by 
phonons can give significant contributions at high tem- 
peratures. In order to have the ballistic regime at suf- 
ficiently low T to minimize the phonon contributions, 
the 2D charge carrier system must have a very high mo- 
bility so that h/ksT becomes very low. Second, the 
temperature constraint, Tp 3> T 3> h/ksT, satisfying 
the dual conditions of being in the ballistic regime (i.e. 
T 3> h/ksT, which is a high-temperature constraint) and 
of also being in the asymptotic low temperature regime of 
T <C Tp (so that the thermal expansion in T /Tp, essen- 
tial for obtaining the linear-T term in the conductivity, 
applies) is not easily satisfied experimentally, and indeed 
most experimentally measured p(T) data in 2D systems 
do not manifest any clear cut linear-T behavior at low 
temperatures. An additional issue we are addressing in 
this work is whether the theory of interaction corrections 
to the conductivity can describe 2D transport in the large 
r s limit as well. This is particularly germane in view of 
the fact that the interaction theory is a systematic many- 
body diagrammatic expansion in the interaction param- 
eter r s (albeit an infinite order formal expansion), and 
the question of the radius of convergence of such an r s - 
expansion becomes quite important for large r s values 
obtained in our samples. 

In this paper, we report our experiments on the low 
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r s near 80 (assuming a hole mass of 0.38m e ). 

In Fig. 0, we show the temperature dependence of the 
resistivity (p) at various densities. For p > 1.7 x 10 10 
cm -2 , p decreases monotonically with decreasing T, 
showing a metallic behavior. For p between 1.2 and 
0.4 x 10 10 cm -2 p shows a nonmonotonic dependence 
on T. It initially increases with decreasing T at high 
T, which-was interpreted as the classical to quantum 
crossovertJ, and then decreases with decreasing T at low T 
showing a metallic behavior as the system goes into the 
degenerate regime. This crossover shifts to lower tem- 
perature with decreasing density and the range where the 
metallic behavior is seen becomes very narrow, especially 
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temperature conductivity and the in-plane magnetore- 
sistancc (MR) of two-dimensional (2D) holes with ex- 
tremely high mobility and very low density [r s ~ 17 to 
80) to study the interaction corrections. From the tem- 
perature dependence of conductivity, we clearly observed 
a temperature region where the the conductivity shows 
a linear dependence on T even in this large r s limit for 
a range of densities in the metallic side of the transition. 
However, the Fermi liquid interaction parameter Ffi , de- 
termined from a comparison of the data with the theory 
by Zala et al.u, shows a surprising non-monotonic depen- 
dence on the carrier density with its value lying between 
—0.5 and —0.7. Fq increases in magnitude with decreas- 
ing density for p > 2 x 10 10 cm" 2 (r s < 22) and then 
decreases with decreasing density for p < 2 x 10 10 cm -2 
(r s > 22), a behavior unexpected from a simple extrap- 
olation of the predicted dependence of Fq on small r s . 
A separate measurement of effective g-factor (g*) from 
the in-plane MR provides a further confirmation of the 
unexpected behavior of Fq . g* decreases with decreasing 
density, consistent with the behavior expected from Fq 
for p < 2 x 10 10 cm" 2 through g* = g b /(l + Fq ) . 

The sample used in this study is a heterojunction 
insulated-gate field-effect transistor (HIGFET) made on 
a (100) surface of GaAsa. A metallic gate, separated by 
an insulator (AlGaAs) from the semiconducting GaAs, is 
used to induce the 2D holes at the interface between the 
GaAs and AlGaAs. Ohmic contacts to the 2D holes are 
made by using a self-aligned contact technique which al- 
lows the diffusion of the contact material under the gate 
region. We would like to emphasize that there is no inten- 
tional doping in the sample and the 2D holes are induced 
by the applied gate voltage. This reduces the scattering 
by ionized impurities so significantly that a very high mo- 
bility can be achieved. The mobility (p) of the sample 
reaches 1.8 x 10 6 cm 2 /Vs at a density (p) of 3.2 x 10 10 
cm~ 2 , which is the highest achieved for 2D holes in this 
low density regime. This high mobility makes h/ksT for 
p = 3.2 to 0.7 x 10 10 cm" 2 range from 16 mK to 80 
mK, low enough that the temperature region where the 
metallic behavior is observed indeed corresponds to the 
ballistic regime while phonon contributions are negligi- 
ble. The extremely high mobility also allows us to mea- 
sure the temperature dependence of conductivity down 
to very low densities reaching p = 1.5 x 10 9 cm~ 2 , with lization in 2D| 
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FIG. 1: p vs T for p 
0.23, and 0.15xl0 10 cm 

p measured as a function of p at different T's. The critical 
density p c is marked by an arrow corresponding to a point 
where p is temperature independent. 



for p = 0.5 and 0.4 x 10 1 



For p < 3 x 10 9 cm 



increases monotonically with decreasing T, exhibiting an 
insulating behavior. To identify the critical density, we 
measured p as a function of p at different temperatures 
and the data is shown in the inset. The critical density 
determined from the crossing point, which shows a tem- 
perature independent resistivity, is p c = 3 x 10 9 cm~ 2 . 
This low critical density, the lowest ever observed in 2D 
systems which exhibit the MIT, is consistent with the 
previous observation by Yoon et alB that the critical den- 
sity becomes lower with decreasing disorder in the sys- 
tem. If we use the hole effective mass m* = 0.38m e , this 
critical density corresponds to r s = 57, which is much 
larger than the r s = 37 predicted for the Wigner crystal- 



In Fig. 0(a), we replot the data for p = 3.2 x 10 10 cm~ 2 
to 0.7 x 10 10 cm -2 as a vs T. The data are scaled by <To, 
the value of a extrapolated to T — 0, and offset by 0.1 
for clarity. The metallic behavior is identified by increas- 
ing a with decreasing T at low T for all these densities. 
Clearly, there is a region where a shows a linear depen- 
dence on T as shown by the best fits in the figure with 
solid lines. We note that r s for these densities ranges 
from 17 to 37 and a shows a linear dependence on T for 
such large r™- To compare our results with the theory 
by Zala et ala, we need to consider several points. First, 
as discussed below phonon contributions are negligible 
in the temperature range where the linear dependence is 
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FIG. 2: (a) a vs T for p = 3.2, 2.2, 1.7, 1.2, 0.9, and 0.7xl0 10 
cm -2 from the top. The data are scaled by the value extrap- 
olated to T = 0, and offset by 0.1 for clarity. The solid line 
for each curve is a linear fit in the low temperature region of 
the data, (b) Fq vs p. 



observed, which is below 200 mK for p = 3.2 x 10 10 cm" 2 
and becomes lower for lower densities. Second, this linear 
region is in the ballistic regime (ksT > Ti/t) since h/hgT 
calculated from cto ranges from 16 mK for p = 3.2 x 10 10 
cm" 2 to 80 mK for p — 0.7 x 10 10 cm -2 . Finally, this re- 
gion is also much lower than the Fermi temperature (Tp) 
of the system, which is 500 mK for the lowest density 
p = 0.7 x 10 10 cm -2 . Zala et al. have pointed out that 
the regime where Fq can be treated as a momentum inde- 
pendent constant is T <C (1 + Fq) 2 Tp, A self-consistency 
check after we have determined Fq approximately (but 
somewhat weakly) satisfies this condition. All these al- 
low a direct comparison of our data with their theory. 
From their theory, the slope of this linear dependence is 
directly related to the Fermi liquid interaction parameter 
Fq by the relation 
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Using to* = 0.38to 6 (which was obtained from cy- 
clotron resonance of high density Id = 5 x 10 11 cm" 2 ) 
2D holes on (100) surface of GaAsE3), we have obtained 
from our data Fq as a function of p. The result is shown 
in Fig. H (b). The value of Fq lies between —0.5 and 
—0.7 in the density range we measured from 3.2 x 10 10 
cm" 2 to 0.7 x 10 10 cm" 2 . For p > 2 x 10 10 cm" 2 
(r s < 22), the magnitude-,increases with decreasing den- 
sity. Proskuryakov et a/.llil have also found that Fq for 2D 
holes in (311) A GaAs/AlGaAs heterostructure increases 



in magnitude with decreasing p for p = 2 to 8 x 10 10 
cm" 2 with values between —0.3 and —0.45. The change 
of Fq per density is similar in both experiments, while the 
magnitude of Fq in our measurement is much larger. We 
note that F^_ibund from experiments on 2D electrons in 
Si-MOSFETt 2 ] also has much smaller magnitude, ranging 
from —0.14 to —0.5 for electron densities 1 to 40 x 10 11 
cm~ 2 _-|While the value found from p-SiGe by Coleridge 
et alE^i is somewhat comparable to ours, between —0.55 
and —0.65. What is surprising in our experiment, which 
explores the much lower density (larger r s ) regime, is 
that the magnitude of Fq does not increase monotoni- 
cally with decreasing density. When the density is de- 
creased below 2 x 10 10 cm" 2 (r s > 22) the magnitude 
of Fq decreases again. This is opposite to the predicted 
dependence of Fq on r s valid for small r s . To our best 
knowledge, the dependence of Fq on r s when r s is large 
has not been calculated theoretically, and it is not possi- 
ble to compare our result with any theoretical predictions 
at this time. 

Experimentally, however, an additional test for this un- 
expected behavior of Fq can be made from the MR mea- 
surements under an in-plane magnetic field (-By). The 
in-plane MR provides a way to measure the effective g- 
factor, which is directly related to Fq by the relation 
g*/gb = 1/(1 + Fq), where gb is the bare g-factor. Fig- 
ure [| (a) shows the in-plane MR measured in our sam- 
ple for various densities at 65 mK. The MR increases 
as exp(B 2 ) at low By and exp(B\\) at high B\\, consis- 
tent with an earlier observation for 2D holes on (311)A 
GaAsEJ. Similarly strong MR-has also been observed for 
2D electrons in Si-MOSFETta and the 2D electrons in 
GaAsEj. It has been established that this crossover from 
low field to high field dependences corresponds to full spin 
polarization of the carriersEZI and its position allows the 
determination of g*. The crossover field B* , determined 
from the position where the second derivative of the p vs 
Bu curve becomes maximum, is marked by an arrow for 
p = 3.2 x 10 10 cm -2 in Fig. (a), and the dependence 
of B* on p is shown in Fig. pf (b). For p > p c , B* de- 
creases linearly with decreasing p (best fit given by the 
solid line), and saturates in the insulating side of the MIT 
for p < p c . This behaviorja also consistent with earlier 
observation by Yoon et alS-S for the 2D holes on (311)A. 
A different way of determining B* , using the inflection 
point between high and low field dependences, yields a 
result within the error bar of this plot, and produces an 
error of less than 15 % in g* . 

For the metallic side, g* determined from the relation 
2Ep = g*/iBB* (where Ep is the Fermi energy and /j,b is 
the Bohr magneton) is shown in Fig. || (c) as a function 
of p. g* decreases monotonically with decreasing p for 
the density range measured. The solid line in Fig|2| (c) is 
the result when the best linear fit in Fig||| (b) is used. Al- 
though a quantitative comparison between g* in Fig |(c) 
and the value of g* expected from Fq cannot be made 
since the bare g-factor gb is not well known for holes 
in GaAs, we can qualitatively compare their density de- 
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FIG. 3: (a) p vs B\\ at T = 65 mK and p = 3.2 
0.9, 0.7, 0.5, 0.4, 0.27, 0.22, and 0.18xl0 10 cm 
bottom. B* is marked by an arrow for p = 3.2 x 10*" cm 
(b) B* vs p. Solid line is a linear fit to the data for metallic 
side, p > p c - (c) g* determined from B* . Solid line is the 
result obtained by the linear fit in (b) 



pendences. From the density dependence of Fq shown in 
Fig. H (b) and the relation g*/gb = 1/(1 + Fq), we expect 
that g* increases with decreasing p for p>2x 10 10 cm~ 2 
and then decreases with decreasing p for p < 2 X 10 10 
cm~ 2 . The behavior of g* shown in Fig ^ (c), therefore, 
agrees well with that expected from Fq for p < 2 x 10 10 
cm -2 , where Fq decreases in magnitude with decreas- 
ing p. While this region of p is where we have observed 
the unexpected density dependence of Fq, in this region 
there is good agreement in the behavior of Fq and g*. 
Thus, our in-plane MR measurements confirm the unex- 
pected behavior of Fq found from the T dependence of 
a. 

The decrease in magnitude of Fq with decreasing p 
(increasing r s ) is of great interest, and needs further ex- 
amination. We, therefore, have also analyzed our data 
without the assumption of a density independent mass 
to* = 0.38TO e . In this analysis, we used three inde- 



pendent relations (equation (1), 2Ep — g*fisB*, and 
g*/gb = 1/(1 + Fq )) between Fq, m* , and g* to calcu- 
late each quantity. Since gb is not well known, we used it 
as a parameter, and found that the density dependence 
of each quantity does not depend on a specific value of gb 
while gb = 0.5 gives the best quantitative agreement with 
Fq and g* determined earlier. From this analysis, Fq 
also exhibits a nonmonotonic dependence on p. The mag- 
nitude of Fq increases with decreasing p for p > 1.5x 10 10 
cm~ 2 and decreases with decreasing p for p < 1.5 x 10 10 
cm~ 2 . The decrease of Fq in magnitude with decreasing 
p for large r s is still observed and confirmed once again. 

Any explanation for this surprising result should take 
into account the large r s values in our system. For 
r s > 37, the 2D system is expected to be a pinned Wigner 
crystal. The critical density for MIT in our system cor- 
responds to r s = 57, considerably larger than the critical 
r s predicted for this crystallization. The r s values for 
which we observed the unexpected behavior of Fq range 
between 22 and 37, close to that predicted for Wigner 
crystallization. Ws note that there has been a Monte 
Carlo calculationta of Fermi-liquid parameters for r s up 
to 5, where Fq still increases monotonically in magni- 
tude with increasing r s . To our knowledge, there is no 
theoretical calculation of the dependence of Fq on r s 
when r s is large, relevant to our experiment. The ques- 
tion whether the crystallization is preceded by a ferro- 
magnetic instability with Fq = — 1 has to be addressed 
as well. -From Monte Carlo calculations, Tanatar and 
CeperleyO have showed both possibilities of a diverging 
and a finite-valued spin susceptibility as r s increases to- 
ward the critical r s for the Wigner crystallization. Our 
result appears to imply that the ferromagnetic instabil- 
ity does not occur in the large r s regime of our 2D hole 
system. 

We now address the important issue of phonon scatter- 
ing contribution to our measured hole resistivity, which 
we have ignored in our analysis. The question of phonon 
contribution to the resistivity is crucial since, if it is sig- 
nificant, it would then be impossible to compare our mea- 
sured resistivity to the interaction theory. We have there- 
fore theoretically directly calculated the phonon scatter- 
ing contribution by including both deformation poten- 
tial and piezoelectric coupling of the 2D holes to GaAs 
acoustic phonons. Following ref. [l9| we have carried out 
a detailed calculation of the phonon scattering contri- 
bution to the hole resistivity in the parameter range of 
our experiment with the results being shown in Fig. ^|. 
Our theoretical phonon-only resistivity, as shown in Fig. 
|], demonstrates that for T < 200mK 1 the temperature 
regime we concentrate on in comparing our experimental 
resistivity with the interaction theory (see Fig. |^), the 
phonon contribution to the resistivity is miniscule (less 
than 1% of the measured resistivity). We are therefore 
justified in neglecting phonon scattering effects in the dis- 
cussion of our experimental results as long as we restrict 
ourselves to T < 200mK as we have done in analyzing 
our data. As is obvious from our theoretical results pre- 
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FIG. 4: Resistivity due to phonon scatterings (piezoelectric 
and deformation potential coupling) as a function of temper- 
ature. Here the lines corresponds to the hole density p — 0.5, 
0.7, 0.9, 1.2, 1.7, 2.2, 3.2xl0 10 cm~ 2 (from top to bottom). 

sented in Fig. |], hole-phonon scattering contribution to 
hole resistivity becomes non-negligible for T > 200mK 
and is, in fact, significant for T > 500mK with its quanti- 
tative importance increasing with decreasing carrier den- 
sity. 

We should point out in this context that we dis- 
agree with the methodology employed recently by 
Proskuryakov et alS3 in subtracting out a calculated 
phonon scattering contribution to their measured resis- 
tivity in analyzing their hole transport data in the con- 
text of a quantitative, comparison with the interaction 
theory of Zala et alu. First, it is well-knownEJ that 
Matthiessen's rule does not apply to 2D systems at fi- 
nite temperatures, and therefore, subtraction of phonon 
contribution (even if this contribution were accurately 
known, which is questionable) in order to obtain the non- 
phonon part is unjustified and may be subject to large 
errors. (This problem is worse in the presence of screen- 
ing of hole-phonon interaction, which must be included 
in the theory.) Second, the calculation of phonon scat- 
tering contribution to hole resistivity, following ref. |2l], 
carried out by Proskuryakov et alx-H is rather crude and 
approximate (compared, for example, with our theoret- 
ical calculations shown in Fig. of this paper). We 



note that in ref. [0] the measured hole resistivity (before 
any phonon subtraction) hardly manifests any clear cut 
linear temperature regime, and the subtracted phonon 
contribution is a large fraction of the measured resistiv- 
ity, thereby casting substantial doubt on the accuracy 
of the subtracted resistivity eventually compared with 
the interaction theory. Our analysis in this work avoids 
these serious pitfalls of ref. [ll] by directly considering the 
measured resistivity in the context of interaction theory, 
which we justify by explicitly calculating the phonon con- 
tribution to the hole resistivity in the temperature range 
of our interest and showing it to be negligible so that 
no arbitrary and unjustifiable phonon subtraction is re- 
quired (in contrast to ref. [ll]). 

We calculate the temperature dependence of the hole 
resistivity by considering screened acoustic-phonon scat- 
tering. We include both deformation potential and piezo- 
electric coupling of the 2D holes to 3D acoustic phonons 
of GaAs. Details of the acoustic-phonon scattering the- 
ory are given in Ref. [l9| In this calculation we use the 
parameters corresponding to GaAs: c/ = 5.14xl0 5 cm/s, 
c t = 3.04 x 10 5 cm/s, p = 5.3 g/cm 3 , eh u = 1.2 x 10 7 
eV/cm, and D = —8.0 eV. In the low temperature range 
(T < 0.2K) we find p(T) cx T 5 - T 7 because the de- 
formation potential scattering dominates over piezoelec- 
tric coupling. In the intermediate temperature range 
(0.2K < T < T B g) we have p(T) cx T 3 mostly due to the 
piezoelectric scattering^!!. Above the Bloch-Gruneisen 
temperature, Tbg = 2kFCi/ks ~ 1 — 2K, both scattering 
processes give rise to linear temperature dependence of 
the resistivity, p(T) cx T. The phonon contribution to 
the resistivity shows very weak temperature dependence 
and is negligible when the temperature is substantially 
below Tbg- We emphasize that the phonon contribution 
to the resistivity cannot be linear for T < Tbg, & n d for 
T < 200mK the phonon contribution is negligible. 

In summary, we have measured the temperature de- 
pendence of the metallic conductivity of extremely high 
mobility dilute 2D holes in GaAs in large r s limit. We 
find that the conductivity shows a linear dependence on 
temperature in the ballistic regime. The Fermi liquid in- 
teraction parameter Fq obtained from our data is found 
to exhibit a nonmonotonic dependence on density and 
decrease in magnitude with increasing r s for r s > 22. 

This work is supported by the NSF and the MRSEC 
at Princeton University, and by US-ONR at Maryland. 
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